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Spin relaxation in undoped quasi-spherical CdS quantum dots at zero magnetic fields is investi- 
gated using time- and polarization-resolved transient absorption measurements. Unlike in previous 
studies of these systems, the measured signals were corrected for spin- insensitive contributions to 
the exciton bleaching dynamics, allowing us to determine the pure spin-related exciton dynamics. 
To explain the observed room-temperature spin-relaxation time of several nanoseconds, we propose 
a novel mechanism based on intralevel exciton transitions with the emission of one LO phonon, 
the absorption of another LO phonon, and an electron spin flip, which is driven by the electron- 
hole exchange interaction. The transition rates, calculated in the present work for different sizes of 
quantum dots and temperatures, are in fair agreement with our experimental results. 

PACS numbers: 72.25.Rb, 78.47.+p, 72.25.Fe, 78.67.Hc, 78.55.Et 



Spin relaxation and decoherence in quantum dots 
(QDs) have attracted increasing attention over the 
last few years. This is i.a. due to the long spin- 
decoherence times observed in QDs, which have led to 
suggestions about their possible applications in quantum 
computation 1,2 . The underlying mechanisms of spin re- 
laxation in semiconductor QDs are still being debated 
(see Ref. 3 for a recent review). The term "quantum 
dot" is actually used for a variety of different objects. 
While most device concepts assume highly symmetric 
QDs, real QDs are usually strongly anisotropic. For ex- 
ample, the most often investigated self-assembled QDs 
have a base elongated along the [110] axis 4 , and QDs 
formed by interface fluctuations in narrow quantum wells 
are elongated along the [110] axis^. As a result, the 
anisotropic exchange interaction splits the | ± 1) radia- 
tive exciton doublet, with opposite total angular momen- 
tum projections, into two linearly polarized eigenstates 
\X) = (|1) + |-1)) /V2 and \Y) = (|1) - | - 1»/(<V2). 
On the other hand, chemically synthesized QDs§>£&£ are 
nearly spherical, retaining higher symmetry, and there- 
fore the exciton states | ± 1) are degenerate eigenstates 
of the Hamiltonian. Research on chemically synthesized 
quasi-spherical QDs is motivated also by their possible 
use as building blocks for constructing artificial solids 
of semiconductor QDs in the bottom-up self-assembly 
approach 10 . However, there exist only few reports on 
spin relaxation in quasi-spherical QDs. Measurements 
of differential transmission at zero magnetic field in neu- 
tral QDs revealed rather small circular polarization of 
the signal, from which only qualitative conclusions could 
be drawr&ii. The results retrieved from the decay of 
the Faraday rotation indicate that the spin dynamics is 
considerably slower in quasi-spherical QDs — than in 
anisotropic QD o 12 ' 13 i 14 . However, as the signal of Fara- 
day rotation contains also contributions from the spin- 
insensitive exciton dynamics, these experiments cannot 



be used for the precise determination of the electron spin- 
relaxation time T\ in quasi-spherical QDs. The majority 
of published works about spin relaxation in QDs concen- 
trate on low-temperature measurements, where most of 
the spin-relaxation processes are strongly suppressed 15 . 
In the present paper, we report on spin relaxation (cor- 
rected for a finite carrier life time) for neutral excitons 
in quasi-spherical CdS QDs with the focus on room tem- 
perature, which is most appealing for potential applica- 
tions. We observed two quite distinctive [(sub)picosecond 
and nanosecond] components in the measured dynamics 
of the circular polarization of the differential transmis- 
sion signal, which are attributed in our model to re- 
laxation processes without and with electron spin flip, 
respectively. The dominant mechanism of the dynam- 
ics, observed at the nanosecond time scale, is identified 
as caused by two-LO-phonon intralevel transitions with 
electron spin flip, driven by the electron-hole exchange 
interaction. 

Experiment. The samples used in our experiments were 
CdS QDs in a glass matrix, a well-known model mate- 
rial of mutually isolated quasi-spherical QDs 9 . Here we 
report results for glass filters made by Hoya and Schott, 
which contain QDs with wurtzite lattice, volume filling 
factor of about 0.1% and typical radii 1.9 nm (sample 
GG435), 2.1 nm (sample Y-44), and 2.4 nm (sample 
Y-46). Our findings are quite general for this type of 
material as was verified by our experiments with sam- 
ples made by different manufacturers (Corning, Toshiba). 
The distribution of the QD sizes (as determined from 
TEM analysis) can be described well by the lognormal 
function (e.g., for sample Y-44 we obtained a distribu- 
tion with average radius 2.1 nm and standard deviation 
0.27). The spin dynamics of the carriers were measured 
using the time- and polarization-resolved transient ab- 
sorption technique of ultrafast laser spectroscopy. Spin 
polarization of optically excited carriers was achieved by 
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absorption of circularly polarized pump pulses and the 
carrier dynamics was deduced from the measured depen- 
dence of the probe-pulse transmission on the time de- 
lay between pump and probe pulses. As a measure of 
transmission changes we used the differential transmis- 
sion AT /T = (T E - T ) /T , where T E (T ) is the trans- 
mission with (without) the pump pulse. A Ti-sapphire 
laser with 82 MHz repetition rate was the source of 80 fs 
pulses. The pump power was kept at a low level, where 
the dynamics is independent of the pump intensity (on 
average much less than one electron-hole pair per QD 
and per laser pulse is photoexcited) . The nonequilib- 
rium spin polarization is reflected in the difference be- 
tween the differential transmission dynamics measured 
with probe pulses co- and counter-polarized with respect 
to the pump pulses (called difference signal in the fol- 
lowing). This difference signal can be directly measured 
using a lock-in amplifier and a photo-clastic modulator, 
which periodically modulates the helicity of the pump 
pulses. However, the difference signal decays not only 
due to carrier spin flips but also because of energy re- 
laxation and recombination of carriers. On the other 
hand, the sum of the differential transmission dynamics 
measured with co- and counter-circularly polarized probe 
pulses (sum signal) is sensitive only to carrier recombi- 
nation and energy relaxation. Therefore, the ratio of the 
difference and sum signals, which is called here circular 
polarization of the signal (Pc), can be used to determine 
the spin-relaxation time T\. 



In Fig. 1 we show typical room-temperature difference 
and sum signals as well as the corresponding dynamics 
of the Pc- The origin of the spin- insensitive dynamics 
(sum signal) was already investigated in details*. There- 
fore, we will concentrate only on the spin-sensitive signal 
in the following. However, the rather pronounced decay 
of the sum signal on the investigated time scale [note the 
logarithmic scale in Fig. 1(a)] implies that correcting for 
the finite carrier lifetime is essential for an accurate deter- 
mination of the spin relaxation time. The decay of the 
Pc on the nanosecond timescale can be well described 
with an exponential decay function. As seen from the 
inset in Fig. 1(b), there is also a (sub)picosecond compo- 
nent of the Pc decay beyond the time overlap of pump 
and probe pulses. As implied by our analysis, this short 
decay component is not connected to the electron spin 
relaxation and we treat it in more detail elsewhere^. 
Here we concentrate mainly on the long decay compo- 
nent. The long component of the Pc decay at 300 K 
can be characterized by an exponential decay function 
with time constant t\ ~ 5 ns (t\ = Ti/2, see Ref. 13). 
At lower temperatures the dynamics is slower and also 
the absolute value of the Pc is larger [see Fig. 1(b)]. In 
fact, below 200 K the decay of the Pc is so slow that 
the characteristic decay time cannot be precisely deter- 
mined within the carrier lifetime. The dependences of 
the transition rate T = X/T\, corresponding to the time 
constant T% of the measured Pc decay, on the size of 
the QDs and on the sample temperature are shown in 
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FIG. 1: (color online) . Spin dynamics in quasi-spherical QDs 
with typical radius of 2.1 nm. (a) Sum and difference of dif- 
ferential transmission signals measured at 300 K using co- and 
counter-circularly polarized probe pulses with respect to the 
circular polarization of the pump pulses; the difference signal 
was normalized to the initial value of the sum signal. The 
inset shows the optical absorption spectrum of the sample, 
the arrow indicates the energy, at which the dynamics were 
measured, (b) Dynamics of Pc at 300 K and 10 K (open 
points) . The curve for 300 K is the exponential decay fit with 
time constant 5 ns. The horizontal line for 10 K is a guide to 
the eye. The inset shows the fast initial decay of the Pc, the 
dotted line illustrates where pump and probe pulses overlap 
in time (as measured by two-photon absorption). 



Fig. 2. The spin- relaxation times, inferred from our mea- 
surements at temperatures above 200 K, are considerably 
larger than those in self-assembled QD o 12 ' 13 and in QDs 
formed by interface fluctuations in QW— . Previously re- 
ported measurements of spin decoherence in chemically 
synthesized quasi-spherical QDs also revealed a nanosec- 
ond component in the decay of the Faraday rotation^^ 
but these measurements were not corrected for the carrier 
lifetime, which left open questions about the relative im- 
portance of the spin-sensitive versus the spin-insensitive 
contributions to the measured decays. In our undoped 
samples the absorption of light leads to the generation 
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FIG. 2: (color online). Dependence of the room temperature 
spin-relaxation rate V on the QD radius 7? (a) and temper- 
ature dependence of F in 2.1 nm QDs (b). Points represent 
the values of T = 1/Ti estimated from the measured Pc de- 
cay and curves show the calculated rate of transitions between 
the lowest dipole-active exciton states (15*, lS3/2)_i/2,3/2 an d 
(IS, 15*3/2)1/2, -3/2- The experimental values of F and the 
corresponding error bars were obtained from the processing 
of multiple data sets measured for each temperature. Inset to 
panel (a): size dependence of the exciton energy levels, cal- 
culated with the Luttinger parameters 71 = 1.09 and 72 = 
0.3421 (energies are counted from the band gap of bulk CdS 
and multiplied by R 2 ). Inset to panel (b): Feynman diagram 
for the transition (IS, 1S s/2 )i/2,f z — > (IS 1 , l<S , 3/2)-i/2,-F i • 
The vertices correspond to the electron-hole spin exchange 
(large circle) and the phonon absorption/emission (small cir- 
cles). The total transition amplitude is a sum of contribu- 
tions, described by six diagrams with different time ordering 
of the vertices. 



of neutral excitons. Therefore, the hyperfine interaction 
with nuclear spins, which is believed to be the dominant 
spin-relaxation mechanism for electrons and positively 
charged excitons (trions) in QDs, is ineffective in our case 
(see, e.g., Ref. 18). 

Theoretical model and discussion. In view of the rela- 
tively small spin-orbit splitting in CdS, A so ~ 70 meV^ 9 , 



the six-band mode l 20 ' 21 is applied to describe hole states 
in QDs. A hole, confined to a spherical QD with a cu- 
bic lattice, possesses a definite eigenvalue of its total an- 
gular momentum F, which is the sum of the hole spin 
J (J = 1/2 for split-off hole bands, J = 3/2 for light 
and heavy holes) and the orbital angular momentum of 
the envelope function L. In the hole wavefunction, com- 
ponents with two different values of the orbital angu- 
lar momentum, L and L + 2, are mixed. We label hole 
states as nXp (cf. Ref. 22, 23), where the index X refers 
to the lowest value of the orbital angular momentum in 
a linear combination, which describes the corresponding 
state (X = S,P,D,F,... for L = 0,1,2,3,..., respec- 
tively). The index n = 1,2,3, .. . labels different states 
with the same X and F in order of increasing energy. 
While within the four-band model the lowest energy level 
of a hole in a spherical QD is always IS3/2, the six-band 
model gives for small CdS QDs, with cubic lattice, the 
lowest energy level lP 3 / 2 2 ° l21 . In QDs with a wurtzite 
lattice, the crystal field splits the energy levels nXp into 
sublevels with definite values of |F z | 24 i 25 , where F z is the 
projection of the total angular momentum of a hole on 
the c-axis of the wurtzite lattice. In bulk CdS the cor- 
responding splitting energy is A = 27 meV— . For the 
conduction band, size quantization is relatively strong, 
and it is sufficient to take into account only the lowest 
energy level, IS", when studying the optical properties 
of QDs near the absorption edge. Therefore, the rele- 
vant states of an electron-hole pair are (lS,nXp) Sz ,p z , 
where s z = ±1/2 is the projection of the electron spin 
on the c-axis. It is important to note that the downward 
shift of the (15, 1S 3 / 2 ) s ^,f^ states due to the electron-hole 
Coulomb interaction in a spherical QD 23 is significantly 
larger than that for the (15, ^-Pj,/2)s z ,f z states. As a re- 
sult, for CdS QDs with R > 1.5 nm the lowest exciton 
states turn out to be (15, 15 3 / 2 ) s f [see the inset to 
Fig. 2(a)]. 

Only for those dipole-active exciton states, for which 
the projection of the total angular momentum on the 
c-axis, N z — F z + s z , equals ±1, the optical absorp- 
tion is influenced by the sense of circular polarization 
of the absorbed light. For the lowest states of this 
type, (15, 153/2)1(11/2, ±3/2, the oscillator strength is 3 
times larger than that for the states (15, 15 3 / 2 )±i/2.±i/2> 
which are separated from the lowest exciton level by 
a relatively wide energy spacing (about 23 meV at 
R ~ 2 nm). Therefore, near the absorption edge, the 
states (15, 153/2)^1/2, ±3/2 are significantly more active 
in the optical absorption compared to higher states with 
N z = ±1, so that the dynamics of the Pc is mainly 
determined by the time evolution of the pump-induced 
populations of the exciton states (15, 15 3 / 2 )_i/2,3/2 an d 
(15, 153/2)1/2,-3/2- This population dynamics is strongly 
affected by the interaction of the excitons with LO 
phonons. This interaction we describe with the Hamil- 

tonian H- mt = J2x (taOa + 7a <4) > where a\ and (a\) 
are the creation (annihilation) operators for phonons of 
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the A-th mode, and j\ are the corresponding interac- 
tion amplitudes derived in Rcf. 23. In spherical QDs, 
both bulk-like and interface phonon modes are charac- 
terized by definite values of the quantum numbers l\ and 
toa, which are related, respectively, to the phonon angu- 
lar momentum and its z-projection. For CdS quantum 
dots with R ~ 2 nm, the finite lifetime of the LO modes 
due to their decay into acoustic vibrations causes an LO- 
mode broadening of about 1 meV— . In general, the mis- 
match between the exciton-level energy spacings and the 
LO phonon energy significantly exceeds this LO-mode 
broadening, so that the probabilities of one-LO-phonon 
transitions between the lowest exciton states are negli- 
gible. In contrast, quasi-elastic two-phonon relaxation 
processes, in which one LO phonon is absorbed and an- 
other is emitted, can lead to fast intralevel relaxation 
of F z (without electron spin flip), contributing signifi- 
cantly to the decay of the Pc on the (sub)picosecond 
time scale. For CdS QDs with R ~ 1.5 to 3 nm, the 
room-temperature rates of the two-phonon transitions 
(IS, IS 3/2 ) Sz ,f z — > (^S,IS 3 / 2 ) Sz ,~f z , calculated within 
the present model, lie in the range 0.2 to 10 ps _1 . Those 
values are comparable with the rates of two-phonon tran- 
sitions between the exciton states \X) and \Y), which 
dominate the decay of the linear polarization of photolu- 
minescence in flat asymmetric InAs/GaAs QDs at tem- 
peratures > 100 K (see Ref. 27). 

However, the Pc dynamics, observed at the nanose- 
cond time scale, cannot be explained by the two-phonon 
transitions discussed above. A key element is that the 
electron spin polarization, induced by a circularly po- 
larized pump pulse, is not affected by these transitions, 
which eventually lead to a reduction of the Pc to ap- 
proximately 1 /3 of its initial value, but not to zero. We 
attribute further decay of the Pc to transitions with elec- 
tron spin flip driven by the electron-hole exchange inter- 
action, which for spherical QDs can be written a o 24 i 25 
Vexch = -f £cxcha^ (r e - r h ) (s • J). Here r e and s (r h 
and J) are the vector coordinates and the spin operator 
of the electron (hole), ao is the lattice constant of the QD, 
and e eX ch is the exchange strength constant, which can be 
estimated for CdS as £ C xch = 35 me V 25 ' 28 . Although the 
exchange interaction is relatively weak (its characteristic 
energy does not exceed 1 meV in spherical CdS QDs with 
i? > 1.5 nm), from the present calculations it is found 
that this interaction provides - in combination with the 



exciton-phonon interaction - an effective channel for exci- 
ton spin relaxation processes, which do require electron 
spin flip. We have identified and calculated the corre- 
sponding Feynman diagrams, which give the amplitude 
of the transition (15, 1S 3 / 2 )i/ 2 ,f z ~ * (^-S,lS 3 / 2 )_i/ 2 _ Fs 
[see inset to Fig. 2(b)]. 

In Fig. 2, we plot the calculated transition rate between 
the lowest dipole-active exciton states (IS, 15 , 3 / 2 )-i/2,3/2 
and (IS, 15*3/2)1/2, -3/2- The calculations performed take 
into account transitions, induced by phonon modes with 
l\, ly = 0, • • • 5 (i.e., by s-, p-, d-, /-, g- and /i-phonons 
of a spherical QD - cf. Ref. 23), through the interme- 
diate exciton states (\S,nS 3 / 2 ), (lS,nP 3 / 2 ), (lS,nS 1 / 2 ), 
(lS,nP 1/2 ), (lS,nD 5/2 ), (lS,nP 5/2 ) with n = l,---6. 
For different sizes of QDs and temperatures, our theo- 
retical results for the transition rate with electron spin 
flip are in line with the measured Pc dynamics. The 
observed trend towards a reduction of T with increas- 
ing R finds its explanation in the fast (roughly as R~ 3 ) 
decrease of the efficiency of the electron-hole exchange 
interaction with increasing R. The temperature depen- 
dence is mainly determined by the decreasing number of 
available LO phonons at lower temperatures. Note that 
no fitting parameters are used in the calculations. 

Conclusions. The dynamics of the Pc in differen- 
tial transmission, measured in the vicinity of room tem- 
perature on quasi-spherical CdS quantum dots, is char- 
acterized by two distinct time constants; one on the 
(sub) picosecond scale, one on the nanosecond scale. We 
show that the Pc decay is governed by two different types 
of transitions between neutral exciton states. While the 
"fast" dynamics is determined by processes, where the 
electron spin is conserved for pump-induced excitons, the 
"slow" dynamics is attributed in our model to electron 
spin relaxation. The mechanism of this relaxation sug- 
gested in the present paper, namely, two-LO-phonon in- 
tralevel transitions with electron spin flip, driven by the 
electron-hole exchange interaction, is shown to provide a 
coherent explanation of our experimental results. 
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